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It is generally difficult to experimentally distinguish magnetic multipolar orders in spin systems. 
Recently, it was proposed that the temperature dependence of the nuclear magnetic resonance 
relaxation rate 1/Ti can involve an indirect, but clear signature of the field-induced spin nematic 
or multipolar Tomonaga-Luttinger (TL) liquid phase [Phys. Rev. B 79, 060406(R) (2009)]. In this 
paper, we evaluate accurately the field and temperature dependence of 1/Ti in spin-i frustrated Ji- 
J2 chains combining field-theoretical techniques with numerical data. Our results demonstrate that 
isotherms of 1/Ti as a function of magnetic field also exhibit distinctive non-monotonic behavior in 
spin nematic TL liquid, in contrast with the standard TL liquid in the spfn-| Heisenberg chain. The 
relevance of our results to quasi one-dimensional edge-sharing cuprate magnets, such as LiCuVC>4, 
is discussed. 
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I. INTRODUCTION 

One of the current topics in solid-state magnetism 
is multiple-spin ordering without any single-spin dipolc 
moment. 1,2 Vector and scalar chiral orders and spin mul- 
tipolar orders are typical examples expected to appear in 
real compounds. The emergence of vector and scalar chi- 
ralities accompanies the spontaneous breakdown of par- 
ity or time-reversal symmetries and these two could be 
detected indirectly from parity- or time-reversal-odd ob- 
servables and related quantities (e.g., electric polariza- 
tion in multifcrroics, asymmetric momentum dependence 
of spin structure factors, Hall conductivity, etc). On the 
other hand, spin multipolar orders, for example, spin ne- 
matic and spin triatic orders, in spin-^ magnets can occur 
without breaking spatial symmetry, since they are char- 
acterized as condensation of bound multi-magnons 3 ~ 8 
or spin-triplet resonating valence bond state. 9 The ab- 
sence of both spin long-range order and lattice symmetry 
breaking makes it difficult to find them in experiments. 
Furthermore, no clear experimental proof of the spin mul- 
tipolar orders has ever been reported, because there is no 
established experimental method of probing spin multi- 
polar orders. 

Thanks to recent theoretical studies, it has been grad- 
ually recognized that magnetic multipolar states occur 
in several geometrically frustrated spin-i magnets, espe- 
cially, in low dimensions. 3 ~ 8 ' 10 ~ 14 In this paper, we fo- 
cus on spin- i frustrated chains with ferromagnetic (FM) 
nearest-neighboring exchange J\ < and competing an- 
tifcrromagnctic (AF) next-nearest-neighboring exchange 
J 2 > 0. The Hamiltonian is given by 

Kji-.h = y] j n Sj ■ Sj +n - h s* , (i) 

n=l,2 j j 

where Sj is the spin-i operator on the jth site, and 



H is the applied magnetic field along the z axis. Re- 
cent studies revealed that this model possesses a series 
of field-induced multipolar Tomonaga-Luttinger (TL) liq- 
uid phases in high magnetization regime. 7 > 8 > 10 > n > 15 These 
multipolar TL liquids are interpreted as a hard-core Bose 
gas of multi-magnon bound states, 7,8 where the num- 
ber p of magnons forming a multi-magnon bound state 
changes successively from p = 2 to 4 with varying J2 / 1 J\ \ . 
In the multipolar phase with p = 2, the quadrupo- 
lar (or spin nematic) operator S^S^ +1 exhibits a quasi 
long-range order. Similarly, the octupolar (or triatic) 
operator S^ 1 S^ +1 S^ +2 an d the hexadecapolar operator 
^j^j+i^j+2^j+3 srl0w a Quasi long-range order, respec- 



tively, in the phases with p = 3 and p = 4. In all of 
these phases, the longitudinal spin correlation also de- 
cays algebraically, while the transverse spin correlation is 
short-ranged. These characteristic properties of the mul- 
tipolar TL liquid phases are distinct from those in the 
standard TL liquid (e.g., spin-i Heisenberg chain). Nu- 
merical studies 8,11 of the quadrupolar and octupolar TL 
liquid phases showed that the relevant multipolar corre- 
lations are dominant in the high-magnetization regime, 
whereas the longitudinal spin-density-wave (SDW) cor- 
relation becomes dominant in the lower-magnetization 
regime. We call the high-field multipolar states with 
p = 2 and 3 quadrupolar and octupolar liquids and the 
low-field multipolar states SDW2 and SDW3 liquids, re- 
spectively, following Ref. 8. 

In these multipolar TL liquids, bulk static quantities, 
such as entropy, specific heat, and uniform susceptibility, 
exhibit qualitatively the same behaviors as those of stan- 
dard TL liquids, and hence it is impossible to distinguish 
multipolar and standard TL liquids from them. To de- 
tect direct evidence of their ordering, we need to measure 
multiple-spin order parameter or relevant multiple-spin 
correlation functions, which is generally difficult. In- 
stead, to find any indirect evidence, it is necessary to 
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specify an effective method of detecting any signature 
of spin multipolar ordering through usual experiments. 
Recently, Ref. 15 has proposed an experimental way of 
detecting a signature of multipolar TL liquids, showing 
that the nuclear magnetic resonance (NMR) relaxation 
rate 1/Xi decreases with lowering temperature T in the 
multipolar TL liquids near saturation. This is completely 
different from the behavior of standard TL liquids in one- 
dimensional (ID) AF magnets such as the spin-i Heisen- 
berg chain, in which 1 jT\ always increases with lowering 
T irrespective of the value of applied magnetic field H. 
However, it is generally hard to experimentally approach 
the vicinity of the saturation field, especially, in magnets 
with strong exchange couplings. In fact, the saturation 
field of a Ji- J 2 magnet LiCuV0 4 is about 40 - 50 T, 16,17 
and to perform NMR measurement under such a high 
field is not an easy task. An experimentally detecting 
scheme at low magnetic field is therefore desirable. It 
was also shown in Ref. 15 that the momentum depen- 
dence of dynamical structure factors in the multipolar 
TL liquids shows a distinct difference from that of the 
spin-i AF Heiscnberg chain. 

In this paper, we reexamine the NMR relaxation 
rate 1/T± in the spin quadrupolar (nematic) and SDW2 
TL liquids, including both low and high-magnetic-field 
regimes, and evaluate its field and temperature depen- 
dence quantitatively. To this end, we combine the field- 
theoretical (bosonization) approach with the density- 
matrix renormalization-group (DMRG) method, substi- 
tuting numerical values for nonuniversal parameters in 
analytic results. We compare these results with the re- 
laxation rate 1/Ti of the spin-i AF Heisenberg chain. It 
is found that, in addition to the temperature dependence, 
the field dependence of 1/Ti in the nematic and SDW2 
liquids clearly differs from that in the usual TL liquid 
of the spin-i AF chain. In the J\-Ji chain, the relax- 
ation rate 1/Ti slowly decreases with increase of magnetic 
field (or is almost independent of the field) in the low- 
field SDW2 regime and increases suddenly in the high- 
field nematic regime near saturation, whereas it increases 
monotonically, function of the field, in the usual 

TL liquid of the spin-^ AF chain. The non-monotonic 
field dependence of 1 /T\ is a unique characteristic of the 
multipolar TL liquids in spin-i J1-J2 chains. Further- 
more, if we tune the direction of the field H, we can 
eliminate the dominant contribution to 1 /T± in the mul- 
tipolar phases and thereby make 1/Ti show exponential 
decays at low temperature in both nematic and SDW2 
liquids. This exponential decay is caused by a spin gap 
in the spin transverse excitations, which is also a unique 
property of the multipolar TL liquids. Our prediction 
would be applicable in NMR experiments of quasi one- 
dimensional (ID) edge-sharing cuprate magnets, such 
as LiCuV0 4 (Refs. 16-19), Rb 2 Cu2Mo 3 Oi2 (Ref. 20), 
PbCuS0 4 (OH) 2 (Refs. 21-24), LiCu 2 2 (Refs. 25-28) 
and NaCu2C>2 (Ref. 29), whose magnetic properties are 
expected to be described by the spin-^ J1-J2 model with 
FM Ji and AF J 2 . 



The paper is organized as follows. In Sees. II and 
III, we review briefly the theory of the NMR relaxation 
rate 1/Ti in electron spin systems and effective theo- 
ries for multipolar and usual TL liquids. Particularly, 
in Sec. Ill, relying on an established field-theoretical 
method, we write down the formula of 1/Ti in multipo- 
lar liquid phases of the J1-J2 chain (1) and also that in 
TL liquid of the spin-i AF Heiscnberg chain. Section IV 
contains the main results of the present paper, which are 
quantitative estimates of 1/Ti m both quadrupolar and 
usual TL liquids as a function of magnetic field and tem- 
perature. Section V is devoted to a discussion of some 
relevant factors in real magnets which are neglected in 
the previous sections. We consider a broad temperature 
scale beyond the low-energy effective theory, the relation 
between electron-nuclear spin interaction and direction 
of field H , and the effects of Dzyaloshinsky-Moriya in- 
teraction. Tuning of the field direction is also discussed. 
Finally, in Sec. VI, we summarize our results and their 
relevance to real quasi-lD compounds. 



II. NMR RELAXATION RATE 

Here, we briefly explain the formula of the NMR relax- 
ation rate 1/Ti in electron-spin systems. Provided that 
the NMR relaxation process is mainly caused by inter- 
action between electron spin Sj and nuclear spin 7, the 
standard perturbation theory for the interaction evalu- 
ates the relaxation rate 1/T 4 as follows: 30,31 

1/Ti ex ]Ti|i ± (fc)| 2 [S+-(k,u) + S- + (k,uj)] 

k 

+ |i||(fc)| 2 5-(fc, W ). (2) 

Here, we have assumed that the electron-spin system is 
in one dimension. The electron-spin dynamical structure 
factor at finite temperature T is defined as 

/oo 
dt e™{S?(t)SS(0)) T , (3) 
j 

where (• • ■ )t denotes the thermal average. The frequency 
u in Eq. (2) is a given resonant value of applied oscillat- 
ing field. Since its magnitude is generally much smaller 
than the energy scale of electron systems, we may take a 
limit oj/k B T = f3u -> in Eq. (2). The symbols A ± (k) 
and .A 1 1 (fc) denote Fourier components of hyperfine cou- 
pling constants between electron and nuclear spins, which 
generally stem from the dipole-dipole interaction and an 
SU(2)-invariant exchange interaction Sj ■ I. The longitu- 
dinal component Ay (fc) originates only from the dipole- 
dipole interaction. In quantum spin systems, A\\ (fc) is 
usually the same order as A±(k). 

The spatial range of interactions between electron 
spins and a single nuclear spin is local, that is, at most 
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the order of the lattice spacing a. Therefore, the k de- 
pendence of Aj_ \\(k) could be negligible. Under the as- 
sumption of such a locality, the hyperfine coupling term 
in the Hamiltonian is given by 



H 



hr 



(4) 



where the electron site closest to the nucleus is assumed 
to be j = and A^ v is the real-space hyperfine coupling 
tensor. In this case, 1/Ti can be approximated as 

1/Ti ex l -Al{j = 0) [S+T H + Sr+( U )] 

+4Cj = 0)S£o( W ). (5) 

Here, hyperfine couplings A±(0) and -A||(0) are propor- 
tional to A xx + A yy and (A 2 :z + A yz y^ 2 , respectively. 
The local dynamical structure factor Sj" (w) at finite 
temperature is represented as 

/oo 
dte^(SS(t)S^0))T. (6) 
-oo 

III. MULTIPOLAR AND USUAL TL LIQUIDS 

In this section, we briefly review the effective theory for 
multipolar TL-liquid phases in the J\-Ji spin chain (1) 
and the well-established TL-liquid theory in the spin-i 
AF Hciscnberg chain. Using these theories, we can derive 
analytic forms of dynamical structure factors and l/Ti. 

A. Theory for multipolar liquid phases 

Two analytic approaches have been developed for field- 
induced multipolar phases in the model (1): One is a 
weak-coupling approach 8 ' 10 where \J\\/ J% C 1 is as- 
sumed and Ji is treated as a perturbation for two de- 
coupled AF-J2 chains. The other is a phenomcnologi- 
cal bosonization in the vicinity of the saturation field. ' 8 
The former weak-coupling theory can lead to only the 
quadrupolar (nematic) phase, but it enables us to cal- 
culate various physical quantities in principle. The 
latter approach can treat all of the multipolar phases 
(quadrupolar, octupolar, and hexadecapolar) , while the 
quantities that can be evaluated are somewhat restricted. 
Here, we employ the latter approach, which is sufficient 
to calculate 1/Ti (the former will be used in Sec. V). 

In the J1-J2 chain with ferromagnetic J\ < 0, the ex- 
citation mode which destabilizes the fully polarized state 
at the saturation field H c is a multimagnon bound state 
consisting of p(> 2) magnons. Below the saturation field, 
the soft bound magnons proliferate to form TL liquid 
with dominant transverse multipolar correlation. The 
system is thus described as a hard-core Bose gas of p- 
magnon bound states. From a numerical calculation of 
the p-magnon excitation energies, it has been shown that 



the soft mode to realize the multipolar TL liquid is two-, 
three-, and four-magnon bound states with momentum 
k = 7T for -2.669 < J 1 /J 2 < 0, -3.514 < Ji/J 2 < 
-2.720, and -3.764 < Ji/J 2 < -3.514, respectively. 7 ' 32 

From the picture of the hard-core Bose gas, we may 
map multipolar operators and the z component of spin 
operators to operators of bosons in the following forms 



J j+P-1' 



bib 



j"3 



(7) 



where bj is an annihilation operator of the hard-core bo- 
son, and we have considered the pth-order multipolar 
liquid states (p = 2, 3, and 4 respectively correspond 
to quadrupolar, octupolar, and hexadecapolar liquids). 
The staggered factor (— 1) J in the first line comes from 
the momentum k — tt of the soft multimagnon bound 
states. To the ID hard-core Bose gas, we can apply the 
standard Abelian bosonization. 33 ' 34 Its low-energy effec- 
tive Hamiltonian is represented as 



H 



eff 

J1-J2 



(8) 



where x = ja is a spatial coordinate, $ and are a 
canonical pair of scalar fields, k is the TL-liquid param- 
eter, and u is the elementary excitation velocity. We set 
K — y 1 here at the saturation limit H — >• H c where the 
hard-core boson gas becomes equal to a free fermion gas. 
It is numerically shown 8,11 that the TL-liquid parameter 
k almost decreases monotonically with lowering magnetic 
field H in quadrupolar and octupolar liquid phases. The 
velocity u is expected to be of the order of | J^i,2 1 except 
for the saturation limit H —> H c where u — > 0. We note 
that the TL-liquid theory is valid only in the temperature 
regime of ksT <C u/a. 

The boson operator bj and the density operator rij are 
expressed in terms of $ and 0: 33 ' 34 



fa + fi cos(V4tt$ + 2npj) + 
-.d x <& + g cos(v / 47r$ + 27rpj) 



(9) 



where p = (1/2 — M)/p is the mean value of the bo- 
son density and M = (SJ) is the magnetization per site. 
Factors /„ and g n are nonuniversal constants. In the 
present representation, the scaling dimensions of e 9lv ^* 
and e q% ^® are, respectively, q 2 n and (? 2 /(4k). 

From Eqs. (7)-(9), two-point correlation functions of 

Sj and p-th multipolar operator = n«=i Sj+ n -i 

at zero temperature are shown to asymptotically behave 
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like 



(S*( T )S*(0)) «M 2 -|-J 



f oofl(2,rpi) 
2 V io 



<WM< p)t (0))o « (-l)Vo 2 (p) 1/(2K) + • • ■ (10) 



in the p-th multipolar TL-liquid phase, where (• • • )o de- 
notes the expectation value of observables at zero tem- 
perature. (Incidentally, the p = 1 case corresponds to 
the usual TL liquid in spin-i AF chain. See the follow- 
ing subsection.) Here, r = it is imaginary time, and 
light-cone coordinates w and w are defined as (iu, to) = 
(x + iur,x — iut). Equation (10) has been numerically 
confirmed, 8,10 ' 11 and this asymptotic behavior survives 
down to the phase boundary between multipolar and 
lower-field vector chiral phases. The power-law behav- 
ior also indicates that in the high-field regime (narrowly- 
defined multipolar regime) with k > 1/2, the multipolar 
correlation is stronger than the longitudinal spin correla- 
tion, while the latter is dominant in the lower-field SDW 
regime with k < 1/2. Namely, k = 1/2 is the crossover 
point between multipolar and SDW regimes. 

In the multipolar TL liquid phase, the transverse 
spin correlation decays exponentially, since a finite en- 
ergy is necessary to violate p-magnon bound states. 
Though we cannot compute the transverse spin corre- 
lation function and the multipolar correlation functions 

(Mj p ' ) (r)M (p ' )t (0))o with p' < p within the theory of 
hard-core Bose gas of multimagnon bound states, we can 
explicitly derive the form of the transverse spin corre- 
lation function (5 l j t (r)S'J(0))o, which is short ranged, 15 
using the weak-coupling theory in the quadrupolar-liquid 
phase (p = 2). Furthermore, the exponential decay of 
the p' (< p) magnon correlations including the trans- 
verse spin correlation has been confirmed numerically 
for the entire regime of the quadrupolar and octupolar 
TL liquid phases. 8 The single-magnon gap can be de- 
termined by measuring transverse dynamical structure 
factors 15 iS ±=F (fc,u;), for example, in neutron scattering 
experiment. It would be at most the order of |Ji| or J 2 . 
This fact of the short-ranged transverse spin correlation 
in the multipolar TL liquids plays a key role in the fol- 
lowing argument. 

The slowly decaying mode in the longitudinal spin cor- 
relator in Eq. (10) contributes to the leading temperature 
dependence of NMR relaxation rate 1/Ti. Through an 
established method based on conformal field theory, 34-38 
we obtain the local dynamical structure factors at finite 
temperature from correlation functions at zero tempera- 



ture and achieve the formula for 1/Ti 



1/Ti cx ^(O^QV 1 



+^f:i(°)y — cos(ktt)B(k, 1 - 2k) 

/27TO 



(11) 



where B(x,y) is the beta function and A» : o(j) = 
= A\\(j). The first and second terms on the 
right-hand side are, respectively, contributions from the 
second and third terms of (S](t)S§(0)) in Eq. (10). If 
we take into account the spatial dependence of A\\ (j) 
instead of the formula (5), A\\. n (j) generally take differ- 
ent values (but they are expected to be the same order). 
The transverse spin correlation also contributes to 1/Ti, 
but it must be a thermal activation type ~ e ~ A /( k BT) ^ 
in which A is proportional to the binding energy of p 
magnons, that is, the single-magnon gap. Therefore, it 
can be negligible when temperature is sufficiently smaller 
than A/ Its- The absence of the contribution from the 
transverse spin correlation is a characteristic feature of 
the multipolar TL liquids. 

The temperature dependence of 1/Ti given in Eq. (11) 
is consistent with the previous prediction in Ref. 15. The 
NMR relaxation rate 1 jT\ diverges with lowering temper- 
ature in the SDW regime (n < 1/2), whereas it decays 
algebraically with lowering T in the high-field multipo- 
lar regime (k > 1/2). Thus the temperature dependence 
changes significantly at the crossover point n = 1/2. This 
feature contrasts sharply with the result of usual TL liq- 
uids (see the next subsection). More detailed properties 
of 1/Ti will be discussed in Sec. IV. 



B. TL-liquid theory for the spin-i AF chain 

Here, we quickly summarize the well-known effective 
theory for the spin-i AF Heisenberg chain in magnetic 
field to compare with the previously mentioned theory 
for multipolar liquids. The Hamiltonian is given by 



(12) 



where J > is AF exchange coupling. As is well known, 
the low-energy physics of this model belongs to the TL- 
liquid universality class from zero field to the saturation 
field H c = 2 J. The low-energy Hamiltonian has the same 
free-boson form as Eq. (8): 



U? = 



Hi 



K-\d x <t>) 2 + K(d x 



(13) 



where 4> and 9 are the pair of scalar fields. The TL-liquid 
parameter K and the velocity v are exactly determined 
from the Bethe ansatz. 39-42 In this model, K changes 
monotonically from 1/2 to 1 when H increases from zero 
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to H c . Similarly to the multipolar TL liquids, the value of 
v is of order of Ja except for the saturation limit H — > H c 
where v — > 0. The TL-liquid theory is valid only for 
k B T <C v/a. 

Spin operators are bosonized as 



Si 



M + —d x (j) + (-l) J ai cos(V47T(?!) + 2irMj) + 



b (-iy + bi cos(V47r<?(> + 2nMj) + 



(14) 



Here ai, bo, and b\ are non-universal coefficients, whose 
values are numerically determined as a function of H (or 
M) in Ref. 43. Equations (13) and (14) enable us to 
estimate 1/Ti. The result is 

K /a\ 2 



i /Tl « 4 (o)-(^) r 1 



(0) y — cos(K%)B(K, 1 - 2i<T) 



27ra 



2if-l 



+^i..(««T"(s) B (s' 

1/(2*0-1 



1 

2~fr 



/27ra 



(15) 



where Aptf) = 4|| : i(j) = A\\{j) and = A ± (j). 

The first and second terms on the right-hand side orig- 
inate from the longitudinal spin correlation, while the 
third term originates from the transverse spin correla- 
tion. In contrast to the multipolar liquids, both longi- 
tudinal and transverse spin correlations yield power-law- 
type functions of temperature. Equation (15) shows that 
cither A» : i or A± : i term diverges at low temperatures 
regardless of the value of K, except for K = 1/2. 44 This 
diverging behavior of 1/T\ with decreasing temperature 
is a common property of usual TL liquids in ID magnets, 
and is quite different from the decaying behavior of 1 jT\ 
in the multipolar liquid phases. 



IV. FIELD AND TEMPERATURE 
DEPENDENCE OF 1/Ti 

Utilizing the formulas (11) and (15), we examine the 
magnetic- field and temperature dependence of 1/Ti in 
both multipolar and usual TL liquid phases. Among 
the three multipolar phases, we concentrate on the 
quadrupolar TL liquid, which is most widely expanded in 
the Ji/ J-2-H space and is believed to be realized in several 
real magnets. We show that 1/Ti in the quadrupolar TL 
liquid exhibits characteristic features distinct from that 
in the usual TL liquid. In Eqs. (11) and (15), the val- 
ues of hyperfine constants A\\j_. n (j) depend on crystal 
structure, the spatial relationship between electron and 
nuclear spins, and the direction of applied field H. For 



simplicity, we assume that all the constants A\\ x-. n (j) are 
equal throughout this section. 

Let us first consider 1/Ti of the quadrupolar liquid 
phase. We can evaluate the expression of 1/Ti given in 
Eq. (11), once the parameters re, u, and go are numeri- 
cally estimated. To this end, we perform numerical calcu- 
lations of spin-i T1-J2 chains using the DMRG method. 
For re and go, we compute the equal-time longitudinal 
spin correlation function (Sj(0)S*i(Q))o and the spin po- 
larization (Sj(0))o in the ground state. We then fit the 
data to an analytic form obtained by the bosonization 
method to estimate re and go- (The details of the method 
have been presented in Rcfs. 8, 43, and 45.) For the ve- 
locity u, we use uniform magnetic susceptibility. From 
the bosonization formula (9), the susceptibility is given 
by 



Xu = d(S*)/dH=p 2 K a/(Tru) 



(16) 



in the p-th order multipolar phase. We can thus obtain 
the value of u from the magnetization curve (i.e., suscep- 
tibility) determined by the DMRG calculation. We de- 
termine the values of re, u, and go at M = 0.05, 0.1, 0.15, 
• • • , 0.45 for several values of J\/J%- We have found that 
the TL-liquid parameter k is basically an increasing func- 
tion of the field H (sec Figs. 14 and 18 of Ref. 8), while 
go has a fairly small field dependence. As expected, u/a 
is shown to be of order of the ID couplings J\ and J2 for 
intermediate magnetization and approaches zero near the 
saturation. For example, in the case of J1/J2 = —1 (—2), 
u/a ~ 1.04J 2 (0.49J 2 ) at M = 0.25, while u/a ~ 0.25J 2 
(O.IJ2) at M = 0.4. 

Figure 1 shows the magnetic-field dependence of 1/Ti 
in the J1-J2 spin chain for J1/J2 = —1 and —2 at sev- 
eral values of temperature. Except for the vicinity of 
the saturation, 1/Ti decreases with increasing H (or M) 
for sufficiently low temperatures fcgT <C u/a ~ \Ji,2\- 
This is mainly due to the monotonic field dependence of 
the TL-liquid parameter k, and is a distinct feature of 
the multipolar TL liquids, in contrast with that of the 
spin-i AF chain (see the next paragraph). This feature 
becomes less significant as temperature increases, and 
it becomes almost independent of the field at relatively 
high temperatures fc^T ~ 0.IJ2. On the other hand, in 
a narrow region near the saturation field, 1/Ti increases 
rapidly with H since the velocity u approaches zero. We 
also confirm the previous prediction 15 on the temperature 
dependence of 1 /T\ : 1 jT\ increases divergently with low- 
ering T in the low-field SDW regime (k < 1/2), whereas 
it decreases algebraically in the high-field quadrupolar 
regime (re > 1/2). Explicit temperature dependence is 
presented in Fig. 2. The characteristic behavior in the 
low-temperature regime ksT < O.IJ2 is consistent with 
the previous statement. 

Next, we compare these features of the relaxation rate 
in the nematic liquid with those in a usual TL liquid of 
the spin- i AF Hcisenberg chain (12). Similarly to the ne- 
matic case, we can evaluate 1/Ti from the formula (15) if 
values for the parameters K, v, a±, and bo are prepared. 
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FIG. 1: (color online) Field dependence of NMR relaxation 
rate 1/Ti in the magnetic quadrupolar TL- liquid phase in the 
J1-J2 spin chain (1) for (a) J1/J2 = —1.0 and (b) J1/J2 = 
—2.0. Left panels show the results for the whole range of 
field < H < H c while right panels show the same data 
in the vicinity of the saturation field in an enlarged scale. 
Solid symbols represent the data obtained from Eq. (11) with 
numerical values of «, u, and go- Lines connecting them are 
guide for the eyes. The vertical dashed lines represent the 
saturation field and the crossover line between the SDW (k < 
1/2) and nematic (k > 1/2) regimes. In (b), the boundary 
to the low-field vector-chiral phase is shown by vertical solid 
line. 



As mentioned in Sec. IIIB, K and v are exactly deter- 
mined by the Bethe ansatz 39-42 and the coefficients a\ 
and 60 have been evaluated numerically. 43 Using these 
accurate values, we numerically determine 1/Ti of the 
model (12). The result is shown in Fig. 3. It indicates 
that 1/Ti is a monotonically increasing (decreasing) func- 
tion of H (T) in the TL-liquid phase of the model (12) 
for sufficiently low temperatures. These two properties 
are completely different from those of the nematic TL- 
liquid phase shown in Figs. 1 and 2. The relaxation rate 
of real quasi- ID spin-i AF magnets has been observed 
in somewhat restricted regime of T and H , for example, 
in Refs. 46-49. Their results seem to agree nicely with 
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FIG. 2: (color online) Temperature dependence of NMR re- 
laxation rate 1/Ti in the magnetic quadrupolar TL-liquid 
phase in the J1-J2 spin chain (1) for (a) J1/J2 = —1.0 and 
(b) J1/J2 = —2.0 and several fixed values of magnetization 
M. 



our results in Fig. 3. 

Finally, let us separately estimate each term in the ex- 
pressions of 1/Ti, Eqs. (11) and (15). They are depicted 
in Fig. 4. Wc find that in the Heisenberg chain (12) the 
relaxation rate 1/Ti is clearly governed by the contri- 
bution from the transverse spin correlation, that is, the 
j4x : i term. In the nematic case, in which the transverse 
spin fluctuation is gapped, the A\\.\ term becomes domi- 
nant in a wide magnetic-field range and is responsible for 
the behavior of 1/Ti decreasing with field H. This fig- 
ure definitely demonstrates that characteristic features of 
1/Ti in the J1-J2 spin chain originates from the absence 
of gapless modes in the transverse spin excitations. 

From these results, we conclude that the NMR relax- 
ation rates 1/Ti in the nematic and standard TL-liquid 
phases exhibit quite different behavior in T-H space, and 
hence we propose that its T and H dependence can be 
used to distinguish the nematic liquid phase from the 
standard TL liquid. Similar characteristic behavior of 
1/Ti is also expected in the higher-order multipolar liq- 
uid phases in the J1-J2 spin chain, since these phases 
share essentially the same properties of spin correlations 
and the field H dependence of the TL-liquid parame- 
ter. Furthermore, a nematic TL-liquid phase in the SDW 
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FIG. 3: (color online) NMR relaxation rate 1/Ti for the stan- 
dard TL- liquid phase in the spin-i Heisenberg chain (12); (a) 
Field dependences for fixed temperatures and (b) temperature 
dependences for fixed magnetization. 
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FIG. 4: (color online) Field dependence of NMR relaxation 
rate 1/Ti for (a) the magnetic quadrupolar TL-liquid phase 
in the J\-3i spin chain (1) with J1/J2 = —1.0 and (b) the 
standard TL-liquid phase in the spin-i Heisenberg chain (12). 
Left (right) panels show the data for k B T I Ji — 0.01 and 
k B T/J = 0.01 (k B T/J 2 = 0.2 and k B T/J = 0.2). In addition 
to the total value of 1/Ti, each term in Eqs. (11) and (15) is 
shown separately. 



regime (k > 1/2) is also shown to appear in a wide pa- 
rameter region of the J1-J2 chain with AF Jx. 45,50 This 
AF-Ji SDW2 phase, where the transverse spin correla- 
tion is short-ranged, is also expected to show a peculiar 
H dependence of 1 jT\ , which is distinct from that in the 
usual TL liquid. 



V. WIDE TEMPERATURE RANGE, 
HYPERFINE COUPLINGS, AND DM 
INTERACTIONS 



So far we have discussed low-temperature behavior in 
spin isotropic pure ID systems. In this section, we con- 
sider the temperature dependence of 1/Ti in wide tem- 
perature range, taking interchain couplings into account. 
We also discuss the relation between the form of hypcrfine 
coupling tensor and the direction of field H, as well 
as effects of Dzyaloshinsky-Moriya (DM) interactions. 



A. Temperature dependence of 1/Ti in a wide 
temperature range 

First we consider a wide temperature window including 
higher and lower temperatures in which simple effective 
Hamiltonians (8) and (13) are no longer valid for describ- 
ing the physics of real quasi-lD magnets. At sufficiently 
low temperature, which is lower than the energy scale 
of weak 3D interchain couplings J3D, long-range-ordered 
phases usually emerge due to 3D couplings and hence ID 
effective theories cannot be applicable. In the case of the 
SDW regime, 1/Ti exhibits a divergence with a critical 
exponent near the critical temperature ~ 1/(T— T c )^, if 
3D ordering occurs through a continuous (second-order) 
transition. In case of nematic regime, the nematic order- 
ing at finite temperature may induce singularity in 1/Ti . 
On the other hand, at high temperature, the effective TL- 
liquid theory becomes unreliable when the energy scale 
kgT is increased up to the order of u/a, which is of the or- 
der of the ID couplings J\ and J2, except for the satura- 
tion limit H — > H c where u — > 0. Deviation from the TL- 
liquid theory also comes from breaking of magnon bound 
states at high temperatures. The binding energy 51 is nu- 
merically estimated, at most, as £"bind — 0.39 J2, 0.28 J2, 
and 0.14J2 for J1/J2 = —2.0,-1.0, and —0.6, respec- 
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FIG. 5: (color online) Schematic behavior of 1/Ti of the ne- 
matic TL liquid in wide temperature region for several pa- 
rameter settings. In the regime Jsd <C ksT <C Jid, our 
prediction based on low-energy effective theories can be ap- 
plicable. 



tively, which gradually decreases with lowering magneti- 
zation and vanishes at the boarder to the lower-field vec- 
tor chiral phase. Thus, our prediction is valid in the tem- 
perature range J^d <C fcsT <C Jid = min[u/a, Sbind]- 
In the higher-temperature regime /c_bT Jid, the re- 
laxation rate converges to a constant value. 52 We draw 
schematic patterns of the temperature dependence of 
1/Ti in the nematic TL liquid in Figs. 5(a) and (b). 



B. Hyperfine coupling tensor and field direction 

Next we discuss the anisotropic case A± ^ A\\ (we 
have assumed Aj_ :Ji = Aii :n so far for simplicity). In the 
extreme case Aj_ ^> An, we can no longer neglect the 
contribution from the exponentially decaying transverse 
spin correlation in 1/Ti. In this extreme case, power- 
law behavior in 1/Ti of multipolar TL liquids is negligi- 
ble down to extremely low temperature, hence we ob- 
serve a thermal activation form ~ e - A /( fe -BT) m \j r £ 1 
in a low-temperature regime fc^T <C Jid- Thus it is 
easy to detect a characteristic feature of multipolar TL 
liquids. If target magnets have sufficiently high crys- 
tallographic symmetry, the principal axes for hyperfine 
coupling tensor A^ v can be defined and the tensor is 



diagonalized. In this case, tuning the field direction par- 
allel to a principal axis, we can eliminate off-diagonal 
elements A xz and A yz , and equivalently set An = 0. 
Therefore, the setup of H parallel to the axis offers an 
easy way of measuring the transverse spin gap and dis- 
tinguishing multipolar liquids from ordinary TL liquid. 
Figures 5(c) and 5(d) present the schematic temperature 
dependence of 1 jT\ in the nematic and SDW2 TL liquids 
for the anisotropic case Aj_ ^> Am . In the opposite limit 
A± <C Am, our predictions in Figs. 1 and 2 are highly 
reliable in J%d <C fcgT <C Jid- However, we should 
note that in this limit Aj_ <C Am, 1/Ti becomes insen- 
sitive to the disappearance of the algebraically-decaying 
transverse spin correlation, the hallmark of the nematic 
liquid, and is not efficient in distinguishing the usual and 
nematic TL liquids. 



C. Effects of DM interaction 

In the rest of this section, we consider the effects 
of magnetic anisotropics which generally exist in real 
magnets. In spin-i systems, one of the most realistic 
anisotropics is the DM interaction defined by 

U D M(Q,<P,m) = ^ COs(Qj + tp)D ■ (Sj X Sj+m) 



(17) 



for m = 1,2. Possible values of Q, <p, and the DM vector 
D depend strongly on the crystal structure of each com- 
pound. Since the DM interaction can sometimes induce 
an excitation gap accompanied with local spin polariza- 
tion, it may violate the multipolar TL liquids. We discuss 
the effects of the DM interaction on the nematic TL liq- 
uid in the low-energy theory. 

In the weak-coupling theory 8,10 ' 53 for the spin-i J1-J2 
chain with | J\ | <C Ji , the nematic phase is expressed by 
the effective Hamiltonian 



-V}cff 

H Jl-J2 



J 7=± 

+d(2a)~ 1 sin(v / 87T0- + ttM) 

+c 2 (d y e+) sin(\/2^6L) + • • • (18) 



under the condition that only the c\ term is relevant 
to the Gaussian model, where y = 2ja, (4>±,0±) = 



±<f> 2 ,0i±02)/V2 > and 



'1(2), » 1(2) < 



is a boson-field 



pair defined in each of the decoupled AF-J2 chains (see 
Sec. IIIB). The coupling constants C1.2 arc propor- 
tional to J\. The ci term makes 0_ pinned. As a 
result, the transverse spin correlation decays in an ex- 
ponential fashion. The remaining (</>+, + ) sector in- 
duces the gapless behavior of longitudinal-spin and ne- 
matic correlation functions. On the other hand, the 
C2 term is known to induce a vector chiral phase 53-55 



with ((Sj x S j+ i) z ) ~ (sin(V27f6L)) ^ in a lower-field 
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regime. 8,53 Let us concentrate on the former nematic TL 
liquid in the following. 

Using the effective theory (18), we investigate five 
typical situations: % DM = Hdm(0,0,1), W}) M = 
H DM (n,0,l), U% M = H DM (0,Q,2), H% m = 
Hdm{tt, 0,2), and W DM = Hdm(tt/2,tt/4,2). Let us 
first consider the case of D = (0, 0, D z ), that is, the DM 
vector is parallel to the applied field H . In this case, 
we find that the bosonized DM couplings contain 
slowly-moving bosonic terms without oscillating factors 



(-iy 



a iq-nMj 



(q is an integer). They are represented 



as 



n DM 



ni in 
n DM 

T/iv 
n DAl 



D z sm(V2^9-) H , (19a) 

Ds{cos(\/27r0-) 2aV2~Trd y 6--\ 

-sin(V2^[(2a) 2 !((W 2 

+ {d y e-f) +•••]}, (19b) 

D zy ^(2a)d y e+ + ■■■ , (19c) 
D z ^{2a)d y 9- + ■■■ . (19d) 



are mvari- 



In general, DM terms with D = (0, 0, D z ) 
ant under a global U(l) spin rotation 5+ 
which corresponds to a shift of the phase field 9 + 
9 + + yj2~pKip. Therefore, the bosonized DM terms with 
D z do not involve any vertex operator with 9 + . In addi- 
tion, all the translationally-symmctric DM terms do not 
contain any vertex term with 4>+ at least for the case of 
incommensurate values of M, since the 2n-site transla- 
tion induces (f>i^(y) —> 4 > i,2(y+2na)+n,y/ir(l/2+M) and 
equivalently <f>+(y) <p+(y + 2na) + nyJVj2(\ + 2M). 
Equation (19) is consistent with these symmetry argu- 
ments. Two DM interactions T-Cj^ M tend to lock 9- 
togcther with the C2 term in Eq. (18), both of which 
compete with the ci term. In the nematic phase, we 
know that the c\ term is most relevant and hence the 
DM interactions I-C^m with sufficiently small SD Z are 
negligible. In the strong DM coupling regime, these DM 
terms compete with the C\ term and could change the ne- 
matic liquid into a chiral phase with ((Sj x Sj + i) z ) ^ 0. 
On the other hand, the leading terms of H^'m can be 
absorbed into the free-boson part by shifting the field 
9±. The shift brings about a finite expectation value 
((Sj x Sj + 2) z ) ~ (dydi^), but we note that these chi- 
ralitics do not accompany any spontaneous symmetry 
breaking. The shift also affects the form of spin dy- 
namical structure factors, that is, the gapless points of 
S ^(k,u>) are slightly changed and a small asymmetry 
of k dependence emerges. In contrast to T~L 1 ^ l l M , do 
not compete with the Ci term. Therefore, chirality and 
nematic quasi long-range order can coexist in the latter 
case of H^'m ■ This symmetry argument indicates that 
the DM terms with D z do not introduce any vertex oper- 
ator in the (0+ , 9+) sector, and this statement would be 
true in the wide nematic-liquid region regardless of the 
value of J\. The gapless nature of the (</>+,#+) sector 



hence survives even after introducing these DM terms. 
From these discussions, we conclude that the nematic 
TL liquid survives even in the presence of DM terms with 
D = (0,0, D z ) at least when D z < |Ji,2|. Thus, the ne- 
matic liquid is relatively stable against DM interactions 
with D = (0,0, D z ). 

For the case of D = (D^O, 0), the U(l) rotational 
symmetry is broken and the emergence of vertex op- 
erators e lq ^ e+ is generally allowed. In fact, W DM is 
bosonized as 



H 



DM 



D x M sin ( y/ir/29 + ) cos ( a 



(20) 



This leading term has scaling dimension 1/(8K + ) + 
1 / (8K- ) and can generate a staggered magnetization 
along the S v axis in each AF-J2 chain. In the weak 
Ji limit, T-L V DM is more relevant than the c\ term with 
scaling dimension 2K- w 2K and it can violate the ne- 
matic liquid phase. [In this limit, K± approaches the 
value of original TL-liquid parameter K (1/2 < K < 1) 
of the AF-J2 chain in magnetic field.] When J\ is 
sufficiently strong, the c\ terms defeats the perturba- 
tion 'H DM and protects the nematic liquid phase, pin- 
ning <p_. In this case, sin(-^/7r/20+) cos(^/ V / '29 '_) gener- 
ates new operators cos(v / 2~7T0+) and cos(v / 2~7T0_) via the 
renormalization-group process. The first term with scal- 
ing dimension 1/(2K + ) can open a gap in the (<f>+, 6+) 
sector and induce a transverse staggered magnetization, 
since the U(l) spin symmetry is broken by the DM term. 
In the weak DM coupling regime, thus, this DM term per- 
turbatively deforms nematic spin liquid, opening a small 
gap in longitudinal spin modes and inducing a small ex- 
pectation value of transverse spins. The other DM in- 
teractions VC^M with D parallel to the x axis do not 
possess any slowly-moving operator within a naive cal- 
culation based on the bosonization, but they generally 
have the ability to generate e lq ^ 9+ , which is allowed 
from the symmetry argument. Thus, the gapless nature 
of the (cp + ,9 + ) sector is expected to be fragile and un- 
stable against DM terms with D ^ (0,0, D z ). 

These discussions on DM terms indicate that if we ap- 
ply magnetic field parallel to the DM vector D, the ne- 
matic TL liquid is stably realized and our prediction of 
1/Ti is reliable in wider H-T space. Even when D is not 
parallel to H, the nematic phase will also survive if the 
DM coupling is small, \Ji^\ S> \D\, and if the system is 
in the temperature regime fcflT> \D\. 



VI. CONCLUSIONS 

We have evaluated accurately the field and tempera- 
ture dependence of the NMR relaxation rate 1 jT\ in mag- 
netic quadrupolar (spin nematic) TL liquid of the spin-i 
J1-J2 chain, combining field-theoretical techniques with 
DMRG results (see Figs. 1 and 2). As a comparison, we 
have also calculated 1/Ti of the spin-i AF Hciscnberg 
chain, using field theories, the Bethe ansatz, and DMRG 
method (see Fig. 3). 



10 



In the nematic and SDW2 TL-liquid phase at a low 
temperature k B T -c l^i.ali the relaxation rate 1/T\ first 
decreases with increasing magnetic field and then rapidly 
increases near saturation. In the higher-temperature 
regime k B T ~ O.IJ2, the field dependence of 1/T\ be- 
comes quite small except in the vicinity of saturation. 
The decreasing behavior of 1 jT\ with increasing H comes 
from the monotonic increase in the TL-liquid parameter 
k, while the rapid increase in 1/T\ near saturation is 
attributed to the decrease of velocity u — > 0. The mono- 
tonic property of k is also essential for the characteristic 
temperature dependence 15 of 1/Ti: with lowering tem- 
perature, l/Xi increases in an algebraic form in the lower- 
field SDW regime (k < 1/2), whereas it decreases in the 
higher-field nematic regime (« > 1/2). These character- 
istic H and T dependences could be a signature of the 
nematic and SDW2 TL-liquid phase. Similar features are 
also expected to appear in higher-order multipolar TL 
liquids, for example, octupolar and hexadecapolar TL 
liquids. Probing the H dependence would be easier than 
doing the T dependence since the former does not nec- 
essarily require the accession to high-field regime, where 
NMR measurements are difficult. A combination of the 
H and T dependences of 1/Ti and the gapless behav- 
ior observed from bulk quantities (such as specific heat 
and magnetic susceptibilities) would present indirect but 
strong evidence for multipolar TL-liquid phases. 

In Sec. V, we also considered the effects of spin 
anisotropics and interchain couplings, which are ne- 
glected in the ideal J\-Ji spin chain model. In particular, 
we point out that if the direction of the external field H 
can be parallel to the principal axis of the hypcrfine- 
coupling tensor, all the algebraic contributions in the 
temperature dependence of 1/Ti, that is, A\\. n terms, dis- 
appear. As a result, 1/Ti of multipolar phases in the J\- 
J2 chain becomes a thermal activation form ~ e ~ A /( k BT) _ 
In addition, we predict that the nematic TL liquid is sta- 
ble for small DM terms with the DM vector parallel to 
the applied field, while it can be easily deformed by the 
DM terms, thereby accompanied with small transverse 
staggered magnetization, when the DM vector is perpen- 
dicular to the field. This indicates that we should apply 
magnetic field parallel to the DM vector to obtain a sta- 
ble nematic liquid phase in real quasi-lD J1-J2 magnets. 

Finally we comment on some real compounds. Re- 
cently, quasi-lD edge-sharing cuprate magnets, for 
example, LiCu 2 02, LiCuVCU, Rb 2 Cu2Mo30i2 and 
PbCuS04(OH)2, have been studied extensively as low- 
dimensional frustrated or multiferroic systems. Their 
magnetic properties are believed to be described by spin- 
2 J1-J2 chains in a certain temperature regime. Ex- 
cept for Rb 2 Cu 2 Mo30i2, a 3D ordering has been ob- 



served below a very low critical temperature T c at least 
in zero magnetic field. Furthermore, the values of J\ 
and J2 have been semi-quantitativcly estimated in several 
ways. The estimated coupling constants and the critical 
temperatures in LiCu 2 2 (Refs. 25-28), PbCuS0 4 (OH) 2 
(Refs. 21-24), and Rb 2 Cu2Mo 3 Oi 2 (Ref. 20) are, re- 
spectively, (Ji/k B ,J 2 /k B ,T c ) ~ (-138 K, 86 K, 24 K), 
- (-13 K, 21 K, 2.8 K), and (Ji/fes, Ji/k B ) ~ (-138 K, 
51 K), while two different results have been reported 
for LiCuV0 4 (Refs. 16-19,56); {Ji/k B , J 2 /k B ,T c ) ~ 
(-19 K, 45 K, 2 K) and - (-182 K, 91 K, 2 K). The 
critical temperatures T c are thus small compared to the 
magnitude of Ji,2 except for LiCu2C>2- It means that 
three dimensionality is small in these compounds. We 
note that our prediction for 1/Ti could be applied to the 
temperature condition k B T c <C k B T <C |^i,2|- 

Among these compounds, the magnetization process 
of LiCuVC>4 has been intensively studied in some exper- 
imental groups. In the low-field regime including zero 
field, spiral phases exist at low temperatures. Above H ss 
7.5 T, the spiral phase turns into another phase which 
was concluded by NMR measurements 19 to be a modu- 
lated collinear phase. The appearance of both the spiral 
and modulated collinear phases was well understood 8 as 
a consequence of the vector chiral phase and the incom- 
mensurate SDW2 phase in the ID J1-J2 spin chain. In 
addition, quite recently, a new phase transition has been 
observed 17 at H sa 40 T, where the saturation field is 
H s w 47 T. Comparing the result with the phase dia- 
gram in the ID J1-J2 spin chain, 8 ' 11 we expect the new 
high-field phase for 40 T < H < 47 T to be a nematic 
long-range ordered phase. So far, the magnetic structure 
of this new phase has not been experimentally identified 
at all. We expect that NMR measurements above the 
critical temperatures of this new high-field phase and of 
the intermediate-field modulated collinear phase would 
be useful to verify whether the new phase is a nematic 
phase. We also note that PbCuS04(OH)2 has a rather 
small saturation field H s « 10 T. 24 This might be an 
ideal material for measuring the H dependence of 1/Ti 
up to the saturation field. 
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